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Abstract 

I address the question of whether supersymmetry provides a viable candidate for the dark matter 
in the Universe. I review the properties of the lightest neutralino as a candidate for solving the 
dark matter problem. I discuss the neutralino's phenomenological and cosmological properties, 
and constraints from present and future experiments. In the minimal supersymmetric model, 
the neutralino mass has been experimentally excluded below some 20 GeV, and is not expected 
to be significantly larger than about 150 GeV. I identify a gaugino-like neutralino as the most 
natural dark matter candidate for a plausible range of parameters. The requirement that the 
lightest neutralino be the dominant matter component in the flat Universe provides non-trivial 
restrictions on other parameters of the model, in particular on the masses of the sfermions. 
Next, I study the consequences of adopting further grand unification assumptions. In both 
scenarios I find sfermion masses most likely beyond the reach of LEP 200 and the Tevatron but 
well within the discovery potential of the SSC and the LHC. I also comment on the effects of 
relaxing grand unification assumptions. Finally, I briefiy outline prospects for the neutralino 
dark matter searches. 



^ Based on invited talks at the XXXII Cracow School on Theoretical Physics, Zakopane, Poland, June 2 - 12, 1992; 
the 23rd Workshop "Properties of SUSY Particles" , Erice, Italy, September 28 - October 4, 1992; and (partly) the 
1993 Aspen Winter Conference on High Energy Physics, January 11 - 16, 1993; to appear in the Proceedings of the 
Erice Workshop. 
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1 Introduction 



Experimentally supersymmetry (SUSY) still remains an open question, many theorists now express 
their growing confidence that it is indeed realized in Nature. Originally applied to high energy 
theories to solve the naturalness problem, supersymmetry has actually proven very successful in 
several other aspects of particle physics ^. Indeed, even minimal supersymmetric grand unified 
model predicts the right ratio of gauge couplings |^ , as recently confirmed by LEP Q ; provides a 
mechanism for dynamical electroweak gauge symmetry breaking |^ ; generically predicts the proton 
decay beyond the present experimental reach j^; and provides a nice connection with theories 
valid at the Planck mass scale, like supergravity and superstrings. At the same time, present 
experimental constraints on even minimal version of supersymmetry, both via direct accelerator 
searches and through supersymmetric loop contributions, are so far relatively mild [Q. It has been 
also recognized for some time 1^ that supersymmetry also provides a nice solution to the dark matter 
(DM) hypothesis which states that perhaps more than 90% of matter in the Universe is non-shining 
(see below and, e.g., Refs. H, [lC|]). 

In this report, I will review the present status of the lightest neutralino, assumed to be the 
lightest supersymmetric particle (LSP), as a DM candidate. The literature on the subject is already 
very extensive and I will have to make some crude selections. I will first summarize the main results 
in the framework of minimal supersymmetry. Next I will emphasize more recent work that has been 
done in the context of the minimal supergravity model. 

In some contrast to more astrophysically and cosmologically oriented reviews on supersymmetric 
dark matter (see, e.g., Ref. [0), in considering the LSP as dark matter I will make a fuller use of 
other, experimental and theoretical, constraints of supersymmetry, including various relationships 
between SUSY parameters. 

I will start by swift reviews of both the DM problem and supersymmetry. In chapter ^ I will 
briefly introduce the minimal supersymmetric model, and will define the lightest supersymmetric 
particle. I will summarize its properties and experimental constraints. In chapter |^ I will discuss 
various cosmological properties of the neutralino and the associated constraints on other supersym- 
metric parameters, in particular on the masses of the sfermions. Next in chapter ^ I will review 
the implications that the requirement that the LSP be the DM in the Universe has on the highly 
constrained grand unification scenario. I will close this chapter with a brief discussion of the im- 
plications of relaxing the grand unification assumptions. In chapter |^ I briefly outline the present 
status of direct and indirect searches for dark matter. Some general comments and a brief summary 
will close this review. 

1.1 The Dark Matter Problem 

There is increasing astrophysical evidence that most of the matter in the Universe is dark, i.e., 
does not emit nor absorb electromagnetic radiation (at least at the detectable level). This evidence, 
while only circumstantial, comes from vastly different cosmological scales, ranging from galactic 
scales of several kiloparsecs (lpc= 3.26 light-year = 3.1 x lO^^m) to clusters of galaxies (several 
megaparsecs), and up to global scales of hundreds of Mpc's. There exist several extensive reviews on 
the subject 0, and I will quote but a few examples. The most well-known evidence comes 

from observations of rotational velocities Wrot(?') of spiral galaxies. If the shining matter were the total 
matter of the galaxies then one would expect Wrot('') l/v^- Instead, one observes Wrot('') ~ const 
which implies that the total mass of a galaxy within radius r, Mdm(?') grows like Mdm(?') oc r. 
One then concludes that there exist extensive galactic halos consisting of DM. Estimates show that 
Mdm/M^is 3 — 10, and perhaps even more. There is somewhat more uncertain evidence on the 
scale of 10 to 50 Mpc (clusters of galaxies) for even more DM. Finally, some first evidence has 
been provided by several groups |l^ (like POTENT ||l^) that over very large scales (hundreds of 
megaparsecs) the total mass density p approaches the critical density pcrit- The critical density 
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Peril = SHq/SttG — 1.9 X 10^^^ {hQ)g/cm^ corresponds to the flat Universe, and ho is the present 
value of the Hubble parameter Hq in units 100 km/s/Mpc. The value fl = p/pcrit = 1 is strongly 
preferred by theory since it is predicted by the models of cosmic inflation and is the only stable 
value for Friedmann-Robertson- Walker cosmologies. Larger values of 57 are also strongly supported 
by present models of primordial barygenesis and most models of large structure formation. 

The visible matter in the Universe accounts for less than 1% of the critical density. Primordial 
nucleosynthesis constrains the allowed range of baryonic matter in the Universe to the range 0.02 < 
ilb < 0.11 ||l^ (and more recently Q fib ~ 0.05). This, along with estimates given above, implies 
that: (i) most baryonic matter in the Universe is invisible to us, and (ii) already in halos of galaxies 
one might need a substantial amount of non-baryonic DM. This is the galactic DM problem. 

If = 1 then most (about 95%) of the matter in the Universe is probably non-baryonic and 
dark. This is the global DM problem. Current estimates show that, for 17 = 1, 0.5 ^ ho ^ 0.7 (the 
upper bound coming from assuming the age of the Universe above 10 bin years), in which case one 
expects the DM abundance in the range 

0.25 <nhl< 0.5. (1) 

More conservatively {i.e., without assuming = 1, but still assuming the age of the Universe above 
10 bin years), one excludes ilhg > 1. On the other hand, the density of DM should be at least 
ilhQ > 0.025 in order to provide minimum required DM at least in galactic halos. 

Finally, one should keep in mind that the uncertainties involved in measuring both Q and ho 
are large and the bounds on flhg quoted above should not be strictly interpreted. (In fact, various 
authors adopt somewhat different criteria - this doesn't cause the conclusions to be dramatically 
different). On the other hand, the evidence that DM is abundant at various length scales in the 
Universe is now so convincing that it would be very surprising if it did not prove true. 

One more remark about the specific nature of DM should be made in closing this section. Before 
the recent COBE discovery of the cosmic microwave background anisotropics iQ (for an elementary 
review, see, e.g., Ref. jl^), models of large structure formation widely favored so-called cold DM, 
rather than hot DM {i.e., non-relativistic and relativistic, respectively, see section |3|). The all-cold- 
DM case, while still very successful, has become increasingly at odds with the measured angular 
correlation function of galaxies and with their pair-wise velocity dispersions. In the aftermath of 
the COBE discovery, some people have argued for a mixed DM scenario with roughly 60-70% of 
cold DM and about 30% of hot DM, like neutrinos with mass in the range of a few eV. (Others 
favor the non-zero cosmological constant as a main contributor to fl = 1.) These issues are far 
from being fully clarified but they won't have any truly dramatic effect on the results presented 
here. Conservatively, all one really assumes by adopting the range (|l|) is that the LSP contributes 
significantly to the total matter density in the flat Universe. 

1.2 Supersymmetry 

Supersymmetry ||l|, ^ has become popular, in part, because it solves one of the most serious problems 
of the Standard Model of electroweak and strong interactions, namely the naturalness and the fine- 
tuning problems, as advocated early, among others, by Veltman [T^ ]. 

Supersymmetry must be broken at some effective scale Msusy since the experimental bounds 
on the masses of the scalar {i.e., spin zero boson) partners of ordinary fermions are experimentally 
known to be higher then the masses of the fermions themselves. The commonly accepted way of 
breaking global supersymmetry is to introduce terms which explicitly break it but in such a way 
that no quadratic divergences are re-generated. Such terms are often called 'soft' and are expected 
to be of the order of Msusy which sets the order of magnitude for the Higgs mass, and therefore 
they should not significantly exceed the Fermi mass scale. Thus in this weakly (or softly) broken 
supersymmetry scenario, the masses of the fermions and gauge bosons are given in terms of the 
Higgs boson's vacuum expectation value, while the masses of the scalars are set roughly by Msusy- 
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In the most commonly accepted approach one starts with a grand unified theory with unbroken 
SUSY above some scale Mqut ~ 10^^ GeV. At that scale, the GUT gauge symmetry breaks down 
to the Standard Model gauge symmetry, but supersymmetry remains unbroken until MguSY ~ 
0(1 TeV). Below that scale, the breakdown of SUSY due to the soft terms induces a non-zero Higgs 
vacuum expectation value, and thus triggers the spontaneous breakdown of the electroweak gauge 
symmetry SU{2)l x U{1)y down to the U{1) of electromagnetism. 

The terms which break global supersymmetry explicitly but softly arise naturally when a su- 
persymmetric GUT is coupled to supergravity. Supergravity is a theory of local supersymmetry 
which provides a framework for incorporating gravity into a unified theory of all four fundamental 
interactions. It also allows for a mechanism of global SUSY breaking. I will come back to this issue 
in section 



4.1 where I discuss the dark matter problem in the context of minimal supergravity. 



The program outlined above is realized even in the simplest phenomenologically viable supersym- 
metric model, known as the minimal supersymmetric model, or Minimal Supersymmetric Standard 
Model. 



2 Minimal Supersymmetric Standard Model 
2.1 Generalities 

The simplest supersymmetric theory of phenomenological (and cosmological) interest is the Minimal 
Supersymmetric Standard Model (MSSM) (for a review, see, e.g., Ref. |], |, 0). In the MSSM, one 
adds to all the fields of the Standard Model their supersymmetric partners to form supermultiplets. 
The supersymmetric part of the Lagrangian results from the following superpotential 

W= {huQLUlH2 + hDQLDlHi + hLLElHi)- fiHiH2. (2) 

generations 

(I use here the same notation for both ordinary fields and their chiral superfields.) Other terms, 
which could break B and/or L number conservation, are absent once a discrete symmetry, the so- 
called i?-parity, is assumed. In Eq. (||) /i(7,d,l are (matrix) Yukawa couplings for the up-type quarks, 
down-type quarks, and charged leptons, respectively, and generation indices have been suppressed. 
The Higgs mass parameter fi is assumed to be of the order of m^. SUSY requires two Higgs doublets: 
Hi = {Hi,H^) and H2 = {H2 ,112). Their neutral components and H2 acquire v.e.v.s vi and 
V2 which give masses to down and up-type fermions, respectively. The Higgs mass spectrum consists 
of two physical scalar fields h and iJ, one pseudoscalar A, and a pair of charged Higgs bosons C^. 
At the tree level the Higgs sector is fully described in terms of just two parameters, which I take 
to be nriA and tan/3 = V2/V1. The expected range of values for tan/3 lies between 1 and rat /rat,. 
Expressions for the Higgs masses can be found, e.g., in Ref. [0. Here I only quote the well-known 
relations: mh < mz\ cos2/3| < mz < mn, rrih < mA < rriH, and mc > mw- Radiative corrections 
to the Higgs masses have recently been shown to be potentially significant. I will discuss their 
implications later. 



2.2 Lightest Supersymmetric Particle 

In the MSSM the four neutralinos Xi* (* = lj---i4) are the physical (mass) superpositions of two 
fermionic partners of the neutral Higgs bosons, called higgsinos and ffj : and of the two neutral 
gauge bosons, called gauginos (bino) and (wino). They are Majorana fermions which means 
that they are invariant under charge conjugation. The neutralino mass matrix is given by j^, 17 
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M2 







Ml 



mz COS 9w COS P —mz sm9w cos f3 
\ —mz cos 9w sin P mz sin 9w sin P 



mz cos 9 w cos P — mz cos 9w sin P \ 
—mzsin9wC0sP mz sin9w sin P 











(3) 



The lightest neutrahno 



X = Xi 



NnW^ + N12B + TViai?? + NuH° 



(4) 



wiU be henceforth assumed to be the hghtest supersymmetric particle (LSP). This assumption is 
a plausible but arbitrary one in the MSSM but is naturally r ealiz ed in the low-energy limit of 
supersymmetric grand unified theories as I will discuss in section 4.1. Due to the assumed i?-parity, 
which assigns a value i? = — 1 to sparticles, and i? = +1 to ordinary particles, the LSP is absolutely 
stable: it cannot decay to anything lighter. It can, however, still annihilate with another sparticle 
(in particular with itself) into ordinary matter. 

The neutralino parameter space is described in terms of four quantities: tan/3 = W2/fi, the 
Higgs/higgsino mass parameter /i, and the two gaugino mass parameters Mi and M2 of the and 
fields, respectively |, |l|, |. 

Usually, one assumes that all gaugino masses are equal at a grand unification theory (GUT) scale 
(see section 4.1), from which it follows that 



Ml = ai/a2M2 ~ O.SMa, 



(5) 



as well as 



M2 = 012I oism— K, 0.3m~ 



(6) 



where ai,2,s are the gauge strengths of the groups C/(l)y, SU{2)l and SU{'i)c, respectively, and m~ 
is the mass of the gluino g, the fermionic partner of the gluon. 

Charginos are the charged counter-partners of neutralinos. Their masses are given by M 0| 



2 1 

m ± ± = — 



Ml +^? + 2m^ 



(7) 



(M| -^J?Y+ Am^ cos2 2/3 + Am^^iMj + + 2M2/isin2/3) , 
and, assuming the relation (^), are described in terms of the same parameters as the neutralinos. 



2.3 Parameter Space of the LSP 

It is first worth recalling the phenomenological structure of the neutralino parameter space. It is 
convenient to display the mass and gaugino/higgsino composition contours in the plane (/i, M2) for 
discrete values of tan/3. This is shown in Fig. 1 for a typical value of tan/3 — 2. For |//| ^ M2, 

w Ml w O.5M2, and the LSP is an almost pure gaugino (and mostly a bino _B"). For M2 ^ 
m^ ~ ImIj Eind the LSP is a nearly pure higgsino. More specifically, it is Hs for fj, < and Hj^ for 
/i > 0, where Hs^a = + ^^2]/\/2- In the intermediate ('mixed') region the LSP consists of 

comparable fractions of both gauginos and higgsinos. To be quantitative, I adopt gaugino purity, 
defined as Pgaugino = -^n + -^12 @, ll^, Hj of, say, 90% to distinguish higgsino, gaugino, and 
'mixed' regions. (For simplicity, below higgsino-like LSPs will often be called just higgsinos, and 
analogously for gauginos.) This distinction is important because cosmological properties of higgsinos 
and gauginos are significantly different! 
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2.4 Present and Future Experimental Constraints on the LSP 

Large fractions of the (/i, M2) plane have been excluded by LEP (see, e.g., Refs. [|l9[ |o|): by direct 
searches for charginos (m ± > 46 GeV) and neutralinos, and indirectly, from their contributions to 
the Z line shape at LEP. If one assumes equal gaugino masses at the GUT scale then additional 
regions are also excluded by the unsuccessful CDF searches for the gluino. An early preliminary 
bound of about 150 GeV (with simplifying assumptions j2l) excludes AI2 ^ 44.7 GeV and leads to 
a lower bound |^ 

> 20 GeV (8) 

on the mass of the lightest neutralino. The bound Eq. (^) also survives possibly large corrections 
to the Higgs masses due to the large top-quark mass. This is because larger statistics accumulated 
during more recent runs of LEP has allowed significantly larger regions of the (/x, M2) plane to be 
excluded also for values of tan f3 as small as one. 

A more careful analysis of the gluino bound which takes into account its possible decays into 
heavier states which next cascade-decay into the LSP leads to a somewhat weaker bound of about 
135 GeV H, or 90 GeV ||, which in turn scales the bound Eq. (|) down to about 18 GeV and 
12 GeV, respectively. 

Recent experimental results also strongly constrain the LSP's properties. First, since the neu- 
tralinos couple to the Z through their higgsino components only, it is not surprising that higgsino-like 
X's below about 45 GeV have been excluded by LEP. Gaugino-like LSPs have only been indirectly 
constrained by chargino and gluino searches, with the result given in the bound (|^). Moreover, the 
regions where the LSP is photino-like (the photino 7 = sinO^W^ + cos O^B is never an exact mass 
eigenstate!) have been excluded above the 96% purity level (l^, |23| as can be seen from fig. 1. 

Significant progress in exploring the neutralino parameter space is expected to be made after 



LEP 200 starts operating in early 1995 26 . The strongest constraint will come from pushing the 
chargino searches up to about 80 GeV, which in turn will place indirect lower bounds of about 40 GeV 
on gaugino-like x's and of about 80 GeV on higgsino-like x's. This will have some cosmological 
consequences |^ as I will discuss later. 

Dark matter searches are also capable of exploring the neutralino parameter space, if one assumes 
the LSP to be dominant component of the galactic halo. I will discuss this in section |^. 



2.5 Theoretical and Cosmological Upper Bounds on the LSP Mass 

While the LSP mass is constrained by high-energy experiments from below, one would also like to 
be able to restrict it from above. 

Theoretically, the expectation that the SUSY breaking scale (and hence also m~) should not 
significantly exceed the 1 TeV scale in order to avoid the hierarchy problem leads, via the GUT 
relation Eq. (||), to a rough upper bound Q ^ 150 GeV, see fig. 1. This upper limit is only 
indicative (and it scales linearly with but sets the overall scale for expected LSP masses. It also 
coincides with a similar bound of about 110 GeV (for mt > 90 GeV) resulting from applying the 
naturalness criterion . (I will have more to say about this in the context of minimal supergravity.) 
This indicative upper bound also implies [l^ that higgsino-like x are strongly disfavored, since they 
correspond to uncomfortably large gluino masses. Experimental constraints and theoretical criteria 
(naturalness) then point us towards the gaugino-like and 'mixed' regions. 

Interestingly, one can derive also cosmological upper bounds on m^. This was originally done by 
Olive, et al. and Griest, et al. in the case of higgsino-like LSP of roughly < 2.8 TeV and 
gaugino-like LSP of about < 550 GeV. These bounds correspond to very large gluino masses of 
about 20 TeV and 3 TeV, respectively. (In the higgsino case the co-annihilation of the LSP with 
charginos and other neutralinos in fact greatly weakens the bound, as I will discuss later.) More 
recently. Drees and Nojiri j2^ have sharpened the bound for the bino case to < 350 GeV in 
the context of minimal supergravity. But they also estimated that, in special cases, (enhanced pole 
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annihilation) in the 'mixed' region the constraint JI/iq < 1 aUows for the LSPs as heavy as 100 TeV 
(!), the range not necessarily fitting the concept of low-energy supersymmetry. One should also 
remember that in deriving such cosmological bounds one assumes the the LSP is a pure higgsino or 
bino state which in practice is never the case. In particular, in the grand-unified scenario, where all 
the masses are inter-related, no well-defined cosmological bound has been found |3^. In brief, it is 
worth remembering that in some cases one can derive cosmological upper bounds on but they 
lie well beyond theoretical expectations. 

3 Neutralino LSP as Dark Matter 

Particle candidates for dark matter can be divided into two broad classes: hot and cold DM pi] . 
The first were relativistic at the time of 'freeze-out' (see below), while the latter were non-relativistic. 
Light massive neutrinos, with the mass of about 10 eV, are the well-known (and well-motivated) 
candidates for hot DM. In the cold DM class the leading candidates are the neutralino LSP Q 
in the mass range of several GeV, and the axion. LEP has excluded several other candidates for 
DM, like heavy neutrinos or their supersymmetric partners, sneutrinos, by placing a lower bound of 
about mz/'2 on their mass. This implies that their contribution to the total mass of the Universe 
is uninterestingly small (fi/ig ^ 10"'^). Other experiments have severely narrowed the allowed mass 
range of the axion In contrast, the lightest neutralino x has all the desired properties of 

a natural candidate for both the lightest supersymmetric particle (LSP) and the dark matter in 
the Universe It is neutral, weakly- interacting, and stable (assuming R-parity). In addition, it 
typically gives the relic abundance in the Universe VlK^ within an order of magnitude from unity. 
In recent years important constraints on the neutralino properties have been derived, and various 
conditions that need to be satisfied for the LSP to be a solution to the DM problem for natural 
ranges of parameters have been pointed out. They will be summarized in this chapter. First I will 
briefly outline the procedure of computing the relic density of neutralinos and, in general, any relic 
species. 

3.1 Computing the Neutralino Relic Abundance 

In the very early Universe {t ^ lO^^^s), all the species, including the neutralinos, were in thermal 
equilibrium with photons [ pO| pT[ . At that stage the neutralinos could reduce their number via 
pair annihilation XX ~^ ordinary matter. As the Universe was expanding and cooling, the density 
of the LSPs was decreasing, and it was becoming increasingly harder for them to find partners to 
annihilate. At some temperature, called the 'freeze-out' temperature T/, the number of the LSPs in 
the Universe became effectively constant, and approximately equal to their number today. In other 
words, the LSP relic density per co-moving volume became approximately constant. 

The relic abundance flhg of the neutralinos, or any other relic species, can be computed by 
solving the Boltzmann (rate) equation 

^ = -3Hn^ {avrei) K - (n^')^] , (9) 

where is the actual number density of the LSPs, and is related to flh^ by flh^ — [m-^n-^hf^) j pcrit] 
n'^ is the number density of the LSP would have in thermal equilibrium at time t; H = R/R is 
the Hubble parameter, R is the cosmic scale factor; a is the cross section for the process xX ~' 
ordinary matter, Vrei is the relative velocity of the two annihilating LSPs, and the symbol () denotes 
a thermal average and sum over spins of the LSPs. It is clear that rate at which decreases is set 
by two factors. One is the expansion of the Universe. The other is due to the excess of the given 
species annihilation over the reverse process. 

Eq. (^ can be integrated either numerically or by means of analytic approximations ||l^, Q . 
It occurs that most methods used in the literature give actually very similar results (up to a few 
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per cent). The scaled freeze-out temperature Xf = Tf/m^ is typically very small (x/ — 0(1/20)), 
justifying at low temperatures the approximation 

{(TVrel) — a + bx (10) 

of the thermally- averaged annihilation cross section (avrei)- The input from high-energy physics is 
thus included in the coefficients a and b. 

Special care must be applied in three cases p^ . This is when the mass of the pair-annihilating 
species is close to: (i) half of the mass of an s-channel-exchanged state (enhanced pole-annihilation); 
(ii) the mass of a new final state (threshold effect); and (in) the mass(es) of other state(s) with 
which a DM candidate can also annihilate (co-annihilation effect). In such cases the expansion ( p^ ) 
fails Q. In particular, the relic density around poles ^ doesn't fall as steeply and is 

larger than a naive estimate would show. The co-annihilation effect has been shown to be of great 
importance for higgsino-like LSPs and will be discussed in section |3.3| . 

3.2 LSP Pair- Annihilation Cross Section 

The DM neutralinos were at the time of 'freeze-out' already very non-relativistic (cold dark matter), 
and therefore they could only pair- annihilate into ordinary matter that was less massive than m^. 

The neutralino relic abundance is thus determined by the annihilation cross section cr(xx 
ordinary matter) of the neutralinos into ordinary matter which in turn depends on many unknown 
parameters: /i. Mi, M2, and tan/3, as well as the top quark, sfcrmions', and Higgs masses. In the 
most general case one needs to specify as many as 27 parameters, including 21 sfermion masses. 
The least number of parameters one has to specify is six: fx, M2, and tan/3 of the pure neutralino 
sector alone, as well as the top mass nit, one of the Higgs masses, e.g., the mass of the pseudoscalar 
niA, and the mass of the sfermions m~, if assumed degenerate for simplicity. With so many free 
parameters to deal with it is thus a priori hard to expect that any solid conclusions can be drawn. 
But in fact this is not the case; that is because, as I will argue later, the LSP relic abundance depends 
primarily on only a few parameters, like fj,, M2 and the mass of the lightest sfermion, and only less 
strongly on the other parameters. (Similarly, in the context of minimal supergravity the number of 



basic parameters is very limited thus allowing for very definite predictions; see section 4.1.) 

Of course, the greater the annihilation rate, the lower the relic abundance. The annihilation 
channels into ordinary fermion pairs, XX ~^ ffi ^re always open, except for the top quark which 
is much heavier. In addition, the following final states may become kinematically allowed at larger 
m^: Higgs-boson pairs, vector-boson Higgs boson pairs, and vector-boson pairs. 

In the region < myy the following annihilation channels can become at some point kine- 
matically allowed. In addition to the mentioned channels //, the LSPs can pair-annihilate into 
hh, AA, hA, hH, AH, Zh and ZA. For > mw new final states open up. They include: W^W~ , 
ZZ, C^C^ , HH, W^C^, ZH, and tt. It occurs that in most of the interesting region of the 
neutralino parameter space, i.e., for both 'mixed' and gaugino LSPs, it is the annihilation into the 
Standard Model fermions that is typically dominant, unless all the sfermion are very heavy. The 
annihilation xX ~^ ff proceeds via the exchange of the Z and the Higgs bosons in the s-channel, 
and via the exchange of the sfermions in the t- and w-channels. Other final states can, however, 
also play an important role and in a full analysis cannot be neglected. The dependence of the relic 
abundance on the parameters involved is significantly different when the LSP is mostly a higgsino, 
or a gaugino, or else a 'mixed' state of both higgsinos and gauginos. This is of course due to dif- 
ferent couplings with which higgsinos and gauginos couple to other matter. A closer study shows 
that the Z exchange is very important in the almost pure higgsino and 'mixed' regions, but not in 
the gaugino region, and the Z pole effect is very broad, unless suppressed by vanishing coupling 
to XX (gaugino and (anti)symmetric higgsino i7(^)5). The exchanges of the Higgs bosons are of 
significantly lesser importance |Q, both because the Higgs bosons couple only to the 'mixed' LSPs 
and because their couplings to the ordinary fermions are suppressed by a factor rrif/mw. Finally, 



9 



the effect of the sfermion exchanges is essentially null for higgsino LSPs, due to a double suppression 
by a factor mf/mw in the invariant amplitude. On the other hand, sfermion exchange in gaugino 
LSP pair-annihilation is dominant, roughly fi/ig oc m~/m^, because the gaugino-fermion-sfermion 

coupling is not suppressed. 

In figs. 2 and 3 I show a few illustrative examples of the LSP relic density in the plane (/i, M2) 
for a typical choice of tan (3 — 2 and niA = 150 GeV. In fig. 2 all sfermion masses are chosen to be 
TOj= 200 GeV. (For > 200 GeV I take it^J = "^xO ^S- «^ ^ show only the case < (the case 
/i > is qualitatively similar) and in the left window all sfermion masses are to— = 400 GeV or , 
whichever is larger. In the right window 1 display the role of one sfermion: 1 keep mj — 400 GeV 
except for one (1 choose the selectron) which is set at to— = 43 GeV which is a current lower limit 
on slepton masses from LEP. These figures will be discussed below. As has been already stated in 
section 1.1, ri/ip > 1 is cosmologically excluded. Assuming the theoretically favored value £1 = 1 
leads to the approximate range 0.25 < ^Ih^ < 0.5. Finally, very small values of relic abundance (like 
fi/ig < 0.025) are probably not cosmologically interesting as they don't solve even the galactic DM 
problem. 



3.3 Higgsinos as DM 

Higgsino- like LSPs with mass below mw, along with gaugino- like ones (but not the 'mixed' states; 
see next section), have been considered excellent candidates for the DM js^, |l^. In the case of 
higgsinos, one could easily find significant regions with the preferred abundance (except near the Z- 
and Higgs poles) for any choice of the sfermion masses mj |l8| as can be seen from figs. 2 

and 3. Above mw, the WW, ZZ, and tt final states cause the relic abundance of higgsino-like LSPs 
to drop well below 0.025, up to a of several hundred GeV [^, |2^ (see figs. 2 and 3 and, e.g., 
fig. 2 of the first paper of Ref. [^). However, as I have mentioned in section from the point of 
view of naturalness higgsinos are not very attractive ||l^ . 

As already mentioned in sect. it has recently been shown that the so-called 'co-annihilation' 
has a devastating effect on the relic density of higgsino-like LSPs, making it uninterestingly small. 
As was originally pointed out by Griest and Seckel , if there exists another mass state x' with 
mass only slightly exceeding the LSP mass, and the cross section for the LSP pair-annihilation is 
suppressed relative to the LSP annihilation with x' then it is this latter annihilation that primarily 
determines the LSP relic abundance. All the conditions for co-annihilation are actually satisfied in 
the higgsino-like region, where typically the next-to-lightest neutralino X2 and the lightest chargino 
xt are only slightly heavier than the LSP. In two recent papers Q (see also Ref. [^) it has 
been estimated that co-annihilation then greatly reduces the higgsino-like LSP relic abundance to 
uninterestingly small values (JI/iq < 0.025). The effect of including co-annihilation on the figures 2 
and 3 is to wipe out all cosmologically interesting domains from the higgsino regions. Thus higgsino- 
like LSPs are neither theoretically nor cosmologically attractive. Clearly, finding a higgsino-like LSP 
at LEP 200, or a chargino that would correspond to the higgsino region, would thus have dramatic 
cosmological implications as it would imply, in the MSSM, a great deficit of DM. 



3.4 'Mixed' LSPs as DM 

Next, 'mixed' LSPs (with comparable fractions of both higgsino and gaugino components) are prob- 
ably not cosmologically attractive either, as they do not provide enough DM in the Universe p^ , 
and often even in galactic halos [Q. Typically flh^ is found below 0.1, or even 0.025. This is true 
both below and above mw JTsf for essentially any choice of independent parameters, like tan /3, mj, 

or Higgs masses. (See figs. 2 and 3.) The smallness of il/ig results from the fact that, in the 'mixed' 
region, some of the LSP annihilation channels always remain unsuppressed up to extremely large 
values of the LSP mass (m^ < 0(100 TeV) pOl). 
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3.5 Gaugino-like LSP: A Natural Candidate for DM 

We are left with the gaugino region 3> M2). Here, the crucially important parameters are the 
sfermion masses. This is because in this region the sfermion exchange in the LSP annihilation into 
// is dominant: II/iq cx m~/m^. This can be seen by comparing figs. 2 and 3. The Z and Higgs 

bosons decouple from gaugino-like neutralinos. For larger also the final states involving Higgs 
boson pairs become important. This happens very roughly when all the sfermions are heavier than 
the heavier Higgs bosons (like H oi A). 

In the MSSM sfermion masses are arbitrary. (For simplicity they are often assumed to be 
degenerate.) However, as was pointed out in Ref. it is the mass of the lightest sfermion (other 
than the sneutrino, whose couplings to gaugino-like LSPs are somewhat weaker, and of the stop 
whose mass is relevant only for > mt) that mostly determines the relic density of the LSPs. In 
other words: the gauginos annihilate most effectively via the lightest charged sfermion exchange. 

The resulting situation can be roughly characterized as follows: 

1. All the sfermion masses are large, 400 GeV. (See the left panel of fig. 3; the graph 
doesn't change much as one varies m~ from 1 TeV down to some 400 GeV.) In this case 

fi/ig > 1 in the region of smaller m^. At larger LSP masses, above the kinematic threshold 
for X annihilation into Higgs-boson final states, the relic abundance drops again, depending on 
the actual masses of the Higgs bosons. (The final state hh, which is always open for > m^, 
is never of any significance.) As one moves from more pure gaugino-like states towards more 
'mixed' states, ilh"^ quickly decreases, and the preferred range 0.25 ^ ri/ip 0-5 forms only 
a relatively narrow strip in the (/i,M2) plane. In the 'mixed' region, Vlh^ remains as always 
below 0.1, and often below 0.025. 

Thus, the all-heavy sfermion scenario is perhaps somewhat less attractive as a solution to the 
DM problem @. 

2. The lightest charged sfermion (say, the selectrone) is neither too heavy nor too light, 120 GeV ^ 
m~< 400 GeV. (Other sfermions could either be also within this mass range, or else heavy. 
Their specific values do, however, influence somewhat the lower bound of 120 GeV: if more 
sfermions are on the lighter side then it goes up since more sfermions contribute to reducing 
fl/ig.) This is cosmologically the most natural range of sfermion masses, for which large do- 
mains of the (/i, M2) parameter space (for gaugino-like LSPs) exist with 0.025 < ilh^ < 1. 
See fig. 2. Within those, there are substantial sub-domains with the preferred abundance 
0.25 ^ HHq < 0.5. Again, the relic abundance decreases towards 'mixed' regions. Also, it 
becomes unacceptably large for asymptotically large values of of several TeV. 

3. The lightest charged sfermion is relatively 'light', 43 GeV ^ rn—^ 120 GeV (and, of course, 
ruy- < m~, "m-j)- (Again, other sleptons could be within this mass range, or else heavier. As for 
squarks, one should keep in mind the current experimental bounds from CDF of 100-150 GeV.) 
See the right panel of fig. 3. In this case, the relic abundance in both the gaugino and 'mixed' 
regions is too small to solve the global DM problem, and often too small (especially for smaller 
values of m~) to even contribute significantly to galactic DM jTs], ^ . If one also includes the 
effects of co-annihilation in the higgsino region, one finds no solution satisfying Eq. (|l]) i n the 



whole (/X, M2) plane. One may thus talk about a cosmological lower bound on m^^ |22], the 
exact value depending somewhat on the masses of the heavier sfermions. 

It is worth stressing that, with co-annihilation included, if even a single slepton were found at 
LEP (below 43 GeVj then the LSP relic abundance would be il/ig ^ 0.07 and thus too small to solve 
the dark matter problem. This can be seen from the right panel of fig. 3. Similarly, if a slepton were 
discovered at LEP 200 then the LSP relic abundance would still be too small to give ~ 1 for 
natural ranges of parameters /i and M2 (roughly below 1 TeV). 
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Finally, one might argue that it doesn't matter if the cosmologically favored region is 'big' or 
'small': even one point in the (/z, M2) plane with large enough relic abundance to solve the DM 
problem (eg. (|^)) should be enough. This is in principle true if the minimal supersymmetry model 
is only considered without its GUT context. However, when the MSSM is viewed as resulting from 
a more fundamental theory valid at the GUT mass scale, several other constraints arise, as I will 
discuss in the next chapter. Thus, even in the MSSM, solutions which require a certain tuning of 
parameters to solve only the DM problem, are in general less attractive, although not excluded. 
This context should be thus kept in mind in judging the cases 1 to 3 as being more or less attractive. 
What I find important, however, is that, despite a large number of free parameters, certain general 
conditions for the LSP to be the dominant matter component in the flat Universe can be established 
and compared with present and future searches for supersymmetric particles at high energy colliders. 



4 The Effects of Grand Unifications 

Until now I have considered the LSP of the minimal supersymmetric model viewed as a phe- 
nomenological theory valid around the Fermi mass scale. The MSSM is clearly motivated by grand 
unifications: one usually allows for only such soft terms that follow from minimal supergravity 
and also assumes that all gauginos share a common (gaugino) mass mi/2 at the unification scale 
(Ml = M2 = rn~ = mi/2) which implies the relations (|^) and (||) at the electroweak scale. These 
assumptions themselves are not a part of the definition of the MSSM (see, e.g., Ref. [0). They 
are, however, usually made in a grand unification scenario, in which the MSSM results from some 
underlying GUT structure valid at the scale Mqut ~ 10"'^^ GeV and is coupled to minimal super- 
gravity. (See, e.g., the review by Ibahez in Ref. Q.) Below I will consider LSP as DM in a fuller 
GUT framework. Next I will comment on the effect of relaxing GUT assumptions. 



4.1 Dark Matter in Minimal Supergravity 

At the GUT scale all the gauge coupUngs become equal |ai = a2 — ctg — acuT ~ Below 
Mqut, the gaugino masses run with the energy scale in the same way as the squares of the gauge 
couplings. In other words, on obtains Mi/ai = M2/a2 = m~/as = 'mi/2/aGUT, from which the 
relations (||) and (^ follow, as well as mi/2 ~ I.2M2. 

In the absence of a commonly accepted 'standard', or 'minimal' grand unified theory, it is usually 
assumed that, at the GUT scale, the MSSM couples to the minimal supergravity model which 
provides the desired form of the supersymmetry soft breaking terms. The minimal supergravity 
model is an effective theory at the GUT scale. Its basic assumptions are described, e.g., in the 
review by Ibanez in Ref. Q (page 56) . I will not review them here but merely quote what is relevant 
to the present report. 

In addition to assuming a unification of gauge couplings and gaugino masses, it is natural to 
expect that at the GUT scale also masses of all the scalars (sfermions and Higgs bosons) are equal 
to some common scalar mass mg. This assumption, in conjunction with SUSY and the gauge 
structure, leads to the following expressions for the masses of the sfermions (except for the top 



squark) at the electroweak scale (see, e.g., Ref. |40 ) 

ml =mi+m?, + br ml,^ ± ml cos2P [T^'''' - Q ^sm"^ 0w], (H) 

Il.r-' jl.r ' 

where /l.r is the left (right) sfermion corresponding to an ordinary left (right) fcrmion, T^""'" and 
Q fh R ^'''2 the third component of the weak isospin and the electric charge of the corresponding 
fermion /, and the coefficients b can be expressed as functions of the gauge couplings at mz and are 
6^6 for quarks and 6 ^ 0.5 for sleptons. 

Thus in the minimal supergravity context there are five fundamental quantities: (a) the common 
gaugino mass 1^1/2, (b) the common scalar mass mg, (c) the Higgs/higgsino mass parameter ^, (d) 
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the common scale Anio of all the trilinear soft SUSY-breaking terms, and (e) the scale BniQ of the 
Higgs soft term (BmofJ-HiH2 + h.c). When the Higgs bosons acquire w.e.w.'s, another independent 
quantity arises which can be chosen to be tan (3 = V2 /vi . The whole spectrum of masses and 
couplings can be parametrized in terms of these six independent quantities. 

It is clear that, in contrast to the MSSM, the masses of the sfermions and the neutralinos are 
now linked to each other via mi/ 2- It is thus interesting to see whether the phenomenological 
and cosmological relations between the neutralino and sfermion masses, previously derived at the 
phenomenological level, can be accommodated within the minimal supergravity context, and whether 
they provide any constraints on the fundamental supergravity parameters mi/2 and mo, and vice 
versa. This issue has been studied from various angles starting with Refs. and more recently 

in Refs. [|3[ ||, |2^, H) and H^. The bottom line is that the requirement that the LSP be the 
dominant DM component in the flat (fl — I) Universe is quite naturally realized in the context of 
minimal supergravity but it also imposes stringent constraints on the fundamental parameters mi/2, 
niQ, and /io (the value of fi at the GUT scale). (More specific results are often sensitive to additional 
assumptions and simplifications made by various authors.) For example, it was found in Ref. | |4^ 
(see also Ref. Q]) that the cosmologically favored region, given by Eq. (P, is only (except near 
poles) realized in the narrow strip of comparable values of mp and mi/2 — I.2M2 (see fig. 4). The 
region toq S> rni/2 was found cosmologically excluded {flhQ > 1). In the region toq ^ '^11/2 there was 
always very little DM {flhQ < 0.025). Cosmologically (i.e., requiring the LSP to be the dominant 
mass component of the fiat Universe) one can derive a lower bound on toq: mo ^ 100 GeV for 
< mvi/citesugradm, and toq ^ 40 GeV in general pc[ |. This, along with other (experimental and 
theoretical) constraints, allows for rather specific predictions at the electroweak energyscale. For 
example, one finds a lower bound on the mass of the lightest slepton of about 100 GeV (see also 
Ref. p6| ), thus confirming a result of Ref. [^ . As for the squarks, which in minimal supergravity 
are invariably heavier, one typically gets to- significantly exceeding 200 GeV |3^, although Drees 
and Nojiri |20 in certain less favored cases (rather large values of Atoq) find m~even below 200 GeV. 

Roberts and Ross have recently performed a careful study of the phenomenological impli- 
cations of the grand unifications with the electroweak gauge symmetry triggered by supersymmetry 
breaking terms. In their analysis, they took into account the corrections to the running of the various 
parameters of the model (gauge couplings and masses) due to multiple mass thresholds above the 
electroweak energy scale. The underlying idea is that above the SUSY breaking scale Msusy one 
uses fully supersymmetric renormalization group equations (RGEs) to evolve the model's parameters 
from the GUT scale down. However, as the masses of the heaviest supersymmetric particles (typ- 
ically the gluino and squarks) become larger than the energy scale Q at which they are evaluated, 
one needs to decouple them from the running of the RGEs. This has to be successively done with 
every new state (mass-threshold), as one further evolves down to the mz scale. These thresholds 
were shown to be of significant importance. Another remarkable feature of their analysis was that 
that they did not fix any of the constraints at the start (like often made to— < 1 TeV) but instead 
compared their relative effects on the final output. On the basis of several independent criteria 
(bounds on and mb, improved naturalness criterion, and improved experimental bounds on as 
and other electroweak and supersymmetric quantities), Roberts and Ross generally found both 
and Too to lie within several TeV, with smaller values (a few hundred GeV) generally favored by the 
values of as and the naturalness criterion. The allowed parameter space, however, still remained 
quite large. 

Their analysis was next extended to include the requirement that the LSP be the dominant 
dark matter component of the flat Universe, see figs. 5 and 6. This constraint was found to be of 
particular importance. First, in the parameter space not experimentally excluded by LEP and the 
Tevatron the lightest neutralino was found to be always the LSP. Furthermore, it invariably came 
out to be gaugino-like, (see fig. 5d) in perfect agreement with previous conclusions ||l8|, 

Remarkably, the dark matter constraint was found to be consistent with other bounds: the 
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favored range of the bottom quark mass (4.15 GeV ^ mi, < 4.35 GeV [Q); the expected range of 
the top mass {implying a rather large value mt Jc; 150 GeV); and the present experimental bounds on 
the strong coupling as{mz) = 0.122 ± 0.010 [Q. All these criteria can be simultaneously satisfied 
(fig. 6) without any excessive fine-tuning of parameters. Furthermore, the dark matter constraint 
eliminated large fractions of the parameter space corresponding to toq 3> (fig- 5c). Also the 

CDF limit rrit > 91 GeV eliminated the region toq <C rni^2 where the LSP is invariably higgsino-like 
for which the relic abundance has been found to be very small. 

The region of the parameter space consistent with all the above constraints is severely limited 
but, remarkably, consistent with supersymmetric masses below 1 TeV. In particular, imposing a 
requirement that the LSP provide most DM in the flat Universe again provides very stringent 
cosmological limits |^ on the masses of the LSP and other particles. One also finds 

60 GeV <m^< 200 GeV, (12) 

the upper limit being also expected in the minimal supersymmetric model ^ on the basis of 
naturalness. Similarly, one obtains |^ 

150 GeV < m+ < 300 GeV (13) 

200 GeV < m~< 500 GeV (14) 

250 GeV < m~ < 850 GeV (15) 

350 GeV <m~< 900 GeV. (16) 

The heavy Higgs bosons are roughly in the mass range between 250 GeV and 700 GeV. Of course, 
lower values of all these masses correspond to less fine tuning and larger values of as ■ The lightest 
Higgs boson tree-level mass invariably comes out close to mz; its one-loop-corrected value Q is 
then roughly in the range 120 to 150 GeV. As a bonus, one also finds 0.116 ^ as{mz) ^ 0.120. 

In summary, the spectrum of supersymmetric particles is expected to lie well within the discovery 
potential of the LHC and the SSC. However, if the LSP is to be the dominant component of matter 
in the flat Universe, then it perhaps less likely, although not impossible, for SUSY to be discovered 
at LEP 200 and/or the Tevatron. 

Of course, the specific numerical results of this and other analyses should be taken with a grain 
of salt. It suffices to say that the findings of recent studies Q are in general consistent 

with each other. More detailed comparisons are, however, almost impossible due to somewhat 
different assumptions and approximations used by various authors. Also, additional effects (like 
radiative corrections to Higgs masses, or mass-thresholds around the GUT scale) and complications 
to the model may alter them somewhat. Nevertheless, it is rather encouraging that the 'minimal' 
set of assumptions along with a number of a priori independent criteria provide self-consistent and 
reasonable predictions for minimal supersymmetry coupled to minimal supergravity. 



4.2 Effect of Relaxing GUT Assumptions 

Finally, one shouldn't forget that the concept of grand unification (and the related assumptions), 
while attractive and in some sense natural, will most likely not be directly tested experimentally 
in the foreseeable future (perhaps never), even though recent LEP results can be interpreted as 
supporting it. Moreover, recent results in superstring theory indicate that the usual GUT assump- 
tions about the common scalar and gaugino masses need not be necessarily hold |5^. One would 
therefore like to learn to what extent the various phenomenological and cosmological results are 
modified when the simple GUT relations given above are relaxed. This question was addressed in 
Rcf. [|l9j (see also Ref. ) in the context of the minimal supersymmetric model without the GUT 
assumptions about the gaugino masses Mi, M2, and m~ (nor of course the supergravity relations 
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between scalar masses). Relaxing GUT assumptions results in potentially significant modifications 
of the phenomenology and cosmology of the neutralino LSP. The study focused on a light neu- 
tralino LSP, within the mass range of a few tens of GeV, since such neutralinos will be searched 
for in accelerators and astrophysical experiments during the next decade. The results can be briefly 
summarized as follows: 

1. As expected, the experimental lower bound ^ 20 GeV doesn't exist anymore, since the 
gluino mass bound from CDF no longer applies; in fact no experimental bound can be derived 
in general. (This is not surprising: it is not always remembered that while experimental 
constraints on charged supersymmetric particles can be regarded as relatively general, the 
corresponding bounds on the neutral particles (like the Higgs bosons or the neutralinos) are 
very model-dependent, or even assumption-dependent. In a sense, experimental lower bounds 
on neutral particles could be regarded as (upper) limits of potential experimental ability to 
search for neutral states in a given model.) However, even though experimentally very light 
neutralinos are now again allowed, it is well known that cosmologically neutralinos below a 
few GeV are forbidden, as they would overdose the Universe (the Zel'dovich bound). 

2. The indicative theoretical upper bound niy^ 150 GeV does not apply either for the same 
reason as above. In principle the LSP could be very heavy, although perhaps one should 
expect rriy- 1 TeV for naturalness reasons. 

3. The composition and mass contours also become significantly modified, depending on how 
different the ratio r = M1/M2 is from the GUT case {rcuT ~ 0.5). Since in the gaugino 
region m-^ « Mi, one finds an interesting scaling property ||l9|| in the (/i, M2) plane: in the 
experimentally allowed region, for smaller values of r (r < tqut), generally lighter neutralinos 
are allowed, whilst for r > rcuT gaugino-like neutralinos are generally heavier than in the 
GUT case. The higgsino region as a whole is also shifted to larger (smaller) values of M2 
for smaller (larger) values of r, but otherwise the mass composition contours are not greatly 
modified. 

4. Gosmologically allowed/favored/excluded regions are also significantly affected by the scaling 
property, depending on the choice of r. For ?- ^ 1, the gaugino-like LSP easily gives the 
closure density. (In contrast, when r ^ 1, the LSP is often predominantly wino-like and, due 
to co-annihilation effects, its relic density is very small [p2[). 

In conclusion, even modest modification of the usual grand unification assumptions has a con- 
siderable impact on the phenomenological and cosmological properties of the LSP. In particular, 
the gaugino-like LSPs as light as a few GeV could still be interesting candidates for DM. Experi- 
mentalists in their search for new physics should not be constrained by theoretical expectations and 
biases. 

5 Direct and Indirect Searches for Supersymmetric Dark 
Matter 

I will only briefly outline the present status and future prospects for the neutralino DM searches. 
More details can be found, e.g., in Ref. |5^. The underlying assumption is that the LSP is a 
dominant component of the halo of our Milky Way, with the relic density of ~ 0.3 GeV/cm^ 
(about 3000 LSPs with mass m^^ = 100 GeV per cubic meter) and the mean velocity of about 
300 km/s. 

One can search for DM neutralinos cither directly, through the halo LSP clastic scattering off 
nuclei ~^ or indirectly, by looking for traces of decays of LSP pair-annihilation products 
(mainly neutrinos). 
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At present direct searches rely on 'warm' germanium and silicon detectors and lack some two order 
of sensitivity necessary to probe the neutralino DM. The problem lies in the fact that neutralinos, 
being Majorana particles, couple to nuclei proportional to their spin (although there also exist 
some coherent interactions with nuclei without spin), and the related scattering cross sections are 
invariably very small. Significant progress is reportedly 54 being made in the R&D of cryogenic 
(T<0.03K) detectors. It is expected that they will achieve desired sensitivity by a dramatic noise 
reduction (mostly /3-radiation and Compton electrons) by measuring the energy deposition (a few 
tens of keV) in both the phonons from nuclear recoil and ionization. First working modules should 
be ready by Spring'93 and a full multi-kg detector by the Fall of 1994. There are also plans to replace 
the presently used spinless ^^Ge with with large-spin nuclei (like ^^Ge with spin=9/2). In Europe, 
another promising and innovative technique based on superconducting detectors is currently being 
developed. (More details can be found, e.g., in Ref. [^ .) 

Indirect searches have been looking for high-energy neutrinos and up-going muons in underground 
detectors. The basic idea can be briefly outlined as follows. Halo LSPs can be captured in celestial 
bodies (Sun, Earth) through their scattering off heavy (core) nuclei. If the scattering is sufficiently 
frequent, they loose enough energy to become gravitationally trapped. After enough LSPs have 
accumulated, they can again start pair-annihilating into leptons, quarks, and gauge and Higgs 
bosons, most of which immediately decay into lighter states. Only neutrinos can make their way 
out of the Sun or Earth core. Their energy is roughly E^, ^ (^ ^ i)'^x- They can pass through 
underground neutrino detectors or, by hitting a nearby rock, produce secondary up-going muons 
that can be searched for in underground muon detectors. Estimates show that for iriy, ^ 80 GeV 
signal from the Sun should be dominant. Lighter neutralinos could more effectively scatter off nuclei 
in the Earth's core (which are of comparable mass) and thus produce a stronger signal. 

It is a bit unfortunate that the LSPs that most effectively interact with nuclei are of the 'mixed' 
type which typically give very little dark matter and thus are not the best candidates for the DM. 
(Their pair-annihilation is in that case also very efficient, see section 3.4.) Nevertheless, some regions 
of the LSP parameter space have already been explored by the 1MB and Frejus groups, and more 
recently by Kamiokande |^5j in searches for up-going muons. Whether or not those regions have 
been excluded depends on what fraction of Galactic halo is shared by the neutralinos. 

Neutrino detectors (like, e.g., MACRO, DUMAND, AMANDA), and a new generation of cryo- 
genic detectors should have a much larger potential to fully cover the neutralino parameter space. 
Much work has, and is being, done both on a technical front and in improving existing theoretical 
calculations of scattering cross sections, capture rates, etc. |Q. First decisive results are expected 
within the next few years. 

A whole spectrum of other methods have been considered. They include looking for products 
of LSP pair-annihilation in the Galactic halo, like monochromatic photons, e"'"-line radiation, and 
in continuum spectrum of cosmic antiproton, positron, and gamma-rays, to mention just of few. 
All these techniques seem to suffer from large theoretical and observational uncertainties and insur- 
mountable backgrounds. At present they are considered less promising. 



6 Agnostic Comments 

Before we all (I mean supersymmetry enthusiasts) become overly excited I would like to remind 
that: 

• Supersymmetry has not been discovered (yet). No 'hints' nor 'evidence', while encouraging, 
cannot be a substitute for a truly unambiguous experimental signal. 

• Dark matter has also not been discovered yet. Same remarks apply. Moreover, even its nature 
(macroscopic objects vs. particles, baryons vs. weakly interacting species, hot vs. cold DM) 
and abundance haven't been definitely established. 



16 



• Cosmic inflation and = 1 still remain (very attractive) hypotheses for which there is at best 
some observational support. 

• The LSP may after all not be completely stable. The i?-parity, invoked to insure the stability 
of the proton, is not a 'fundamental' symmetry, and can easily be relaxed. For example, one 
can allow for lepton or baryon number violating interactions (but not both), in which case the 
LSP would be unstable. Despite stringent limits, such a scenario may be easily realized and 
in fact may be favored in light of the COBE discovery. 

• Simple (simplistic?) assumptions at the GUT scale can easily be modified or relaxed, even 
in the context of superstrings, and thus may lead to significantly altered conclusions at the 
electroweak scale. 

• Minimal SUSY may not be the end of the story! After all, it is not favored by any particular 
reason other than simplicity. On the other hand, it is encouraging that even the simplest 
supersymmetric model is both theoretically and cosmologically attractive, while not running 
into any immediate conflict with experiment. 

• Finally, one may argue that still lacking is a satisfactory explanation of why the LSP relic 
abundance should be close to unity in the first place. (See, however, Ref. ||ll|.) Why/how did 
supersymmetric parameters 'conspire' to make the LSP dominate the Universe? This question 
should be answered (or maybe shown ill-posed?) by a truly fundamental theory. 



7 Final Comments 

In this review I have addressed the issue of whether supersymmetry can provide an attractive 
candidate for solving the dark matter problem. I have argued that the lightest neutralino of the 
minimal supersymmetry has all the desired properties for being both the LSP and DM. This is 
certainly encouraging given the fact that supersymmetry is the leading candidate for the extension 
of the Standard Model. In addition, the experimental values of the gauge coupling, as measured at 
LEP, when evaluated in the minimal supersymmetric model, become equal at the SUSY GUT scale 
(of around 10"'^^ GeV), thus supporting the idea of supersymmetric unification. 

Cosmology provides additional constraints on the supersymmetric parameter space. Requiring 
that the LSP provide enough (or at least a substantial fraction of) DM in the flat Universe, when 
combined with the usual theoretical assumptions and present experimental bounds, points towards 
the gaugino in the mass range of several tens of GeV as the most natural candidate for both the LSP 
and the DM. Significant fractions of this mass range should be accessible to LEP 200 and planned 
experiments for DM searches. In addition, cosmologically favored masses for sleptons and squarks 
lie beyond the reach of LEP 200 and the Tevatron (what may be somewhat discouraging) but well 
within the reach of the future hadronic supercolliders: the LHC and the SSC. During the next 
decade, accelerator and astrophysical DM searches should be able to test this and other predictions 
of supersymmetry. 

It is not unlikely that most of the Universe is actually made of supersymmetric particles while 
what we usually call 'ordinary' matter may look more like a very exotic component of the cosmic 
zoo. (One might look at this as yet another step in the Copernican revolution!) Nothing short 
of a direct SUSY discovery will convince the whole physics community but we should certainly be 
encouraged to continue vigorous searches. It is not unlikely that the first signal of supersymmetry 
may actually come from searches for dark matter. 
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Figure Captions 



Figure 1: Contours of the neutralino mass and gaugino-higgsino compositions (purity, as defined 
in the text) in the (/i, M2) plane for tan/3 = 2 and for /x < (left panel) and /it > (right 
panel). Lightest neutralino mass contours are labeled (in GeV). Contours of constant gaugino 
purity {pgaugino = + Z12) arc shown with Pgaugino = 0.99, 0.9, 0.5, 0.1, and 0.01 from larger \n\ 
to larger M2. (Note Pgaugino — 0.01 implies a 99% higgsino.) In the regions labeled B and limited 
by the dashed curves p~ > 0.99 (almost pure bino), while the region (labeled Hs) above the other 
dashed curve correspond to almost pure symmetric higgsino {phiggsino > 0.99). (The almost pure 
anti-symmetric higgsino Ha lies above the range of M2 in the window /Lt > 0.) The areas marked 
"LEP" are ruled out by LEP and the areas below the curves marked "CDF" are ruled out by CDF. 
Also marked is the remaining photino region with 95% purity. 

Figure 2: Contours of relic abundance in the (/U, M2) plane, with tan/3 = 2, niA = 150 GeV, 
m~ = 200 GeV or (whichever is larger) for all sfermions, and /U < (/i > 0) in the left (right) 
panel. Areas excluded by LEP (CDF) are marked "LEP" ("CDF"). The relic density increases with 
increasing grayness. The cosmologically excluded region Qh^ > 1 is delineated by a thick solid line 
in the grey region. The cosmologically favored regions (0.25 < O/ig ^ 0.5) are delineated by dotted 
lines. The long-dashed and thick-solid contours correspond to 0.1 and 0.025, respectively The LSP 
co-annihilation has not been included. Its effect would be to greatly reduce the relic abundance in 
the higgsino region (M2 |/u|). 

Figure 3: Contours of relic abundance in the (/i,M2) plane, with tan/? = 2, m,A = 150 GeV, 
and /Lt < 0. (The case /x > is qualitatively similar.) In the left panel all nij = 400 GeV or 
(whichever is larger). In the right panel all sfermion masses are as before but for the selectron 
m-— 45 GeV or rn^ (whichever is larger). All the textures are as in fig. 2. The LSP co-annihilation 
has not been included. Its effect would be to greatly reduce the relic abundance in the higgsino 
region (M2 > l/x]). 

Figure 4: Contours of relic LSP density in the (/i,M2) plane for tan/3 = 2, niA = 200 GeV, and 
/i = =f500 GeV in the left (right) panel. The cosmologically favored regions with 0.25 < Qh^ < 0.5 
are shaded, while regions where O/iq > 1 are cross-hatched, and regions excluded by LEP and CDF 
are hatched. The HHq = 0.1 contour is shown as a thin dashed line. Above the lines "1" rrig > 1 
TeV for at least one squark, above "2" > 1 TeV, and to the left of "3" rUq < 100 GeV for at 
least one squark. 

Figure 5: In the plane (^1/2, ?72o) for the fixed ratio fia/mQ = 21 show: in window a) the mass 
contours of the top and the bottom quarks (solid and short-dashed lines, respectively); in window 
b) the contours of as{mz) (solid) and the measure c of fine-tuning (dots); in window c) the relic 
abundance nftp of the LSP; and in window d) the mass contours of the LSP (solid) and the lightest 
chargino (dashed) at 50, 100, 150, 200, 500, and 1000 GeV, starting from left, and the contribution 
(dots) of the bino to the LSP composition (bino purity). In all the windows thick solid lines delineate 
regions experimentally excluded by the CDF (marked CDF) where rrit < 91 GeV and by the LEP 
experiments (LEP) where the lightest chargino is lighter than 46 GeV. In window c) we also mark 
by O/iq > 1 the region cosmologically excluded (too young Universe). The thin band between the 
thick dashed lines in window c) corresponds to the flat Universe (J7 = 1). In window d) the region 
excluded by CDF almost coincides with the bino purity of 50% or less. 

Figure 6: I show a blow-up of the down-left portion of the plane {mi/2,mo) from the previous 

figure for the same fixed ratio jio/mo = 2. I combine the mass contours of the top and the bottom 
quarks with the ones of the LSP relic mass density. I use the same textures as in Fig. 5 but I 
also show (two medium-thick short-dashed lines) the contours mt = 4.15 GeV and 4.35 GeV which 
reflect the currently favoured range of the mass of the bottom quark (see text). I sec that they 
cross the cosmologically favored region (thick long-dashed lines) marked Q = 1 at roughly 150 
GeV < mi/2, mo < 400 GeV and for rrit broadly between 150 GeV and 180 GeV. 
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